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Abstract 

In this paper we firstly obtain two kinds of effective spectral density functions by setting the 
cut-off frequencies of baths be infinite and finite. Secondly, we investigate the reduced dynamics 
of open qubits in four kinds of systems constructed with the basic spin-boson model. It is shown 
that the qubit has different dynamics governed by the two kinds of spectral density functions. 
In addition, we obtained that a qubit coupled to an intermediate harmonic oscillator has longer 
decoherence and relaxation times as they are coupled to a common bath than to their respective 
baths. In solving the dynamics of qubits we use a numerically exact algorithm, iterative tensor 
multiplication algorithm based on the quasiadiabatic propagator path integral scheme. 
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I. INTRODUCTION 



As a quantum system of interest interacts with its environment which is also considered 
being made up of a large number of the quantum particles, solving the dynamics of the total 
system is in fact an impossible task by now. However, in many cases we are only interested 
in the dynamics of the quantum system of interest, thus we can investigate the reduced 
system by using the phenomenological environment as proposed by Caldeira and Leggett 
l) . Usually, in this method the environment is modelled by a bath of harmonic oscillators 
(or other quantum particles) and the influence of the bath on the quantum system is fully 
enclosed in the so-called spectral density function J (00). According to the characteristic of 
the environment one can choose the J(uS) Ohmic, sub-Ohmic, super-Ohmic forms and so 
on. Thus, the reduced dynamics of the quantum systems can be investigated in detail. 
Many such problems have been treated extensively and some significative results have been 
obtained in the field in past years {2I, ^j. 

The spectral density functions J (00) obtained from the microcosmic dynamical model 
plus its phenomenological environment (MDMPPE) is in good agreement with that ob- 
tained from quasi-classical equations of motion satisfied by the system of interest. We can 
obtain the effective spectral density function (ESDF) following Leggett as follows. Setting 

cj5 (lu — ujj) /rriiUi = Jo(u), taking cut-off frequency of the bath modes infinity, and sup- 
posing the spectral density function J(u) a continuous function, we can obtain the ESDF 
with J(u>) =lm[K(u)} (see Ref.j^]), where K(u) represents the effects of arbitrary linear 
dissipative and/or reactive elements. Usually, K(u) can be obtained from the equation of 
motion of the system in Fourier-transformed form, namely, K(u)q (u) = — [dV (q) /dq] (u) , 
where V (q) is a conservative potential. This equation can be obtained from the scheme of 
MDMPPE. The ESDFs from infinite cut-off frequency of baths for many models are obtained 
5KZJ. This method is in fact a universal one. 

However, when the temperature of the systems are low enough, we can imagine that 
the dynamics of the quantum systems of interest obtained from the MDMPPE should have 
some difference with that from quasi-classical equations of motions. In this case, the high- 
frequency modes of the bath are in fact suppressed, so the cut-off frequency may be not 
very large and it may be comparable to the characteristic frequency of the quantum system 
of interest. So there are at least two differences between the dynamics of normal and low 
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temperature for the quantum systems of interest. The first is that the ESDFs obtained from 
the MDMPPE at low temperature should not be in accordance with the spectral density 
functions from quasi-classical equations of motion, and the second is that the system at low 
temperature must go through a non-Markovian process 8|, |9| other than the Markovian one. 

With the progress of technique many quantum systems have been designed and manufac- 
tured. In particular, many qubit models for making quantum computers in future are pro- 
posed in recent years. Most of these quantum systems should work in very low temperature 



. . 1 lH . The non-Markovian effects of baths in the systems have been considered recently 
-15|. It is reported that in many chemical and biologic systems, the non-Markovian ef- 



fects may also be important 



16l-ll9|. At very low temperature, the cut-off frequencies of the 



bath should be finite. Taking a finite cut-off frequency of the bath rather than the infinite 
one, we shall obtain another kind of the spectral density function. In the following we de- 
note the ESDFs from infinite, and from finite cut-off frequencies of the bath for model x by 
Jl((jj) (from infinite cut-off frequencies) and J^(oj) (from finite cut-off frequencies), where x 



denotes the physical models. A somewhat simi 



classical framework has been given in Refs 



20 



ar discussion on the problem in completely 



2l|. 



In this paper, we shall investigate the ESDFs of four models by using the scheme of 
MDMPPE at infinite and finite cut-off frequencies of baths. Through the two sets of ESDFs 
J T X (u) and J x (cj) we shall investigate the reduced dynamics of qubits in some practical 
physical models. They are a quantum system in a bath (model A); a quantum system coupled 
to an intermediate harmonic oscillator (IHO) and only the latter coupled to a bath (model 
B); a quantum system coupled to an IHO and both of them coupled to their independent 
baths (model C) and a common bath (model D). 



II. TWO KINDS OF SPECTRAL DENSITY FUNCTIONS AND THE DYNAMICS 
OF OPEN QUBITS 

In this section, we shall obtain two kinds of ESDFs for the models A, B, C, and D when 
they are reduced to the original spin-boson form. One kind of them are obtained from taking 
finite cut-off frequencies of baths and the other is gotten by setting the cut-off frequencies 
of baths infinity. After obtaining the spectral density functions we shall solve the reduced 
dynamics for the quantum systems of interest when they are set as the qubits. The reduced 
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dynamics governed by the two kinds of ESDFs for the qubits will be compared. 

In order to compare the dynamics of the qubit governed by the different effective spectral 
functions in each model, we choose a same starting point Jo (w) for all of the models. 
In another word, we always reduce them to the original spin-boson model with the same 
IHO transformation and set Jq{uj) = J^i !^ 10 ~ u i) / m i w i = ^we _W//Wc in our derivations. 
Because the model A is just the original spin-boson model, our work for reducing this model 
is just reducing it to itself. So in the assumption of infinite cut-off frequency, the obtained 
ESDF of the model should be equal to the Jq (oj), which will be shown in the following 
subsection A. Models B, C, and D are different from the original spin-boson model, so when 
they are reduced to the original spin-boson form, the ESDFs should be different from the 
starting function J (u). 

We have calculated the dynamics of the quantum system of interest governed by J], (uj) 
and j£ with the same values of r\ in each models of the paper, and they are certainly 
different. However, even for the same model, the different ESDFs (lj) and J% (co), may 
provide different damping for the same r] in the function Jq(uj), which will affect the dynamics 
of the quantum system of interest. Do the different kinds ESDFs affect the dynamics of 
quantum system of interest mainly through resulting in different damping or through their 
different forms? To answer this problem, we introduce an effective coupling strength rf = 



221 ] . The effective coupling constant is 



J(uj = u )/u to calculate the dynamics as Ref. 
introduced from the physical deduce that the modes which has the same frequency as the 
characteristic frequency of the system ujq has the major contribution to the system-bath 
interaction. Thus, by using the effective coupling constant we can obtain that, if the reduced 
dynamics governed by (cu) and j£ [co) are different, we can safely attribute the different 
forms of the J x (u) and (u). For lack of space, in the following, we only plot the evolutions 
of elements of reduced density matrix in the same effective coupling strength rf for each 
models. 



A. Two basic models (models A and B) 

A quantum system embeds in a bath (model A) and a quantum system interacts with 
an IHO which is coupled to a bath (model B), are two basic models and many complicated 
models are developed from them. So in this subsection, we use them to show the Leggett's 
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method for obtaining the ESDFs J\j B (w), and J^/ B (u>). Then from the two sets of ESDFs 
we investigate the reduced dynamics of qubits in the two models. 

Model A: Because the physical quantities of the quantum system of interest is not included 
in the ESDF, we can deduce the ESDFs of the model A with the Hamiltonian as jsj 



2/i ' " ^' ' ~«2mi ' 2"~ 1 " 1 V~' ' rriiufj ' ^ 
where p q is the momentum conjugate to coordinate q. Defining U' = dU/dt and using the 
dots for the time derivatives, we can obtain the classical equations of the motion as 



m = -u'(g) - - E, 



1 ' 

rriiuf 



rriiXi = -rriiUtXi - dq. 



Using the Fourier transforms, we can write Eqs.(j2]) and ([3]) as 

-^u?q (co) = -U'Jq (t)) - J2i°i x i («) - £,^r^> 



—miLo 2 Xi (co) = —vriiLolxi (u) — c^q (u) . 



From Eq.(jSJ) we have 



Xi (u) 



rrii (u 2 — u 2 ) 



q(u) 



Insetting Eq.^BJ) into Eq.flU), we have 



2 



Setting 



and 



we then have 



rriiUjf (uf - u) 2 ) 



q(oo) = -K(q(t)). 



rriiUJ 2 {uj 2 — u 2 )' 



J (J) = vr V -^-8 (u/ - Ui) = wV~ w7wc , 



i poo 

L(u) = — 



J (u/) 



-duo' 



7T J Q CO' (uJ l2 — W) 

1 [°° expHy/a^) ,_ 1 



7T 7q (u/ 2 — W 2 ) 



7T 



where (see Appendix) 
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(2) 
(3) 

(4) 
(5) 

(6) 

(7) 

(8) 
(9) 



(10) 



(11) 



The expression of W (cu) is given in the Appendix. Substituting Eq. CTTj) into Eq.((7]) we have 
K (uu) q (ou) = [-/iw 2 + cuijW (ou) + wr/ze^ /aJC ] q (cu) 

= -K (?) • (12) 

If taking cu c be infinite, we can easily obtain the spectral density function as 

4 (cu) = ur]e-" /wc . (13) 

This is the conventional Ohmic spectral density function of bath in the original spin-boson 
model and it has been widely used in quantum dissipative systems. Here, we have only 
shown that, starting from J (cu), an ESDF J\ (cu) can be obtained. The J (cu) and J* A (cu) 
are the same for the bath in the original spin-boson model. The result shows that the 
Leggett's method give a consistent result for the model. If the model is not the original 
spin-boson model (as following models B, C, and D), it can also be reduced to the original 
spin-boson model but the ESDFs should not be equal to the Jo (^) even we start from the 
same J (cu). 

If taking cu c be finite, we then obtain the ESDF as 

J%(cu) = lm(K(cu)) = cur]& (cu) . (14) 

Here, (cu) = Im (W (cu)) +e~ ul ^ c . The other forms of the spectral density functions, such as 
sub-Ohmic and super-Ohmic forms can be obtained similarly. It shows that when the cut-off 
frequency is finite, we should add the term ?7Co>Ihi[W (cu)} in the spectral density function 
as Eq. ffT4"l) . It is interesting that how does the added term influence the spectral density 
function? In Fig.l we plot Im[W (cu)} in different cut-off frequencies of the bath. It shows 
that the absolute values of the Im[W (cu)] decreases with the increase of cu c . When cu c —> 00, 
the Im[iy (uu)} — > 0. However, when the cut-off frequency of the bath is not big enough, the 
influence of Im[W (cu)} should not be neglected in the description of the environment. 

In above derivation of the ESDFs, the physical quantities for the quantum system of 
interest have not been involved, so the J (cu) can control the dynamics of not only the 
harmonic oscillator as in Eq.((T]) but also other quantum systems. In order to check the 
differences of the influences that above two different ESDFs J\ (cu) and (cu) impose on 
the dynamics for the systems of interest, we firstly investigate a simple quantum system, 
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qubit. For the open qubit the total Harmiltonian reads 

h . pf 1 2 / da 



^-t^ + A^ + E^ + ^^+^j ■ (15) 

This is one of the simplest open models of quantum system but its reduced dynamics can still 
not be solved exactly, and some approximations must be appealed to. Many investigations 



on this topic have been made in recent years 23M27T] . As we stressed that when the cut-off 



frequency is not too large to compare to the characteristic frequency of the qubit, the non- 
Markovian effect is non-negligible. The well established iterative tensor multiplication (ITM) 
algorithm based on the quasiadiabatic propagator path integral (QUAPI) scheme is a good 
tool for solving the reduced dynamics of quantum system of interest. The ITM algorithm 
is a numerically exact one and is successfully tested and adopted in various problems of 



open quantum systems 



works 



28 



2814321] . For details of the scheme, we refer readers to previous 



321 ] . To make the calculations converge we use the time step St = 0.1/ A which 



is shorter than the correlation time of the bath and the characteristic time of the qubit. In 
order to include as much non-Markovian effect of the bath as we can, and avoid the heavy 
calculation load we choose Sk max = 3. Here, Sk max x St is the memory times being included 
in the calculations. In this paper we set T = 300 k, and the initial state of the bath be the 
thermal equilibrium state, namely, pbath(0) = e~ l3Hb /Tr(e~^ Hb ) . Here, (3 = 1/ksT, ks is 
the Boltzman constant. As we calculate the evolutions of diagonal elements of the reduced 
density matrix pn(t), we set A = 1 x 10 12 Hz, e = 0.01A, and the initial state of the qubit 
be pi (0) = |$o) (^ol- As we calculate the evolutions of non-diagonal elements of the reduced 
density matrix puit), we set e = 1 x 10 12 Hz, A = O.Ole, and the initial state of the qubit 
be p 2 (0) = |^o) (^ol- Here, |$ ) = |0) and |tf > = l/v 7 ^ (|0) + |1)). The |0) and |1) are the 
basis states of the qubit. The comparison pictures on pn(t) and p^if) hi different ESDFs 
(J A (u), and J^(cj)) and different cut-off frequencies of the bath {u c = 4. OA, 4.1A, 4.3A, and 
10. OA) are plotted in Fig.2, where 7/ is set to be 0.004. 

It is shown that in the model A the qubit has shorter decoherence and relaxation times 
governed by the Ja(u>) than by J a {oj). 

Model B: There is another basic model that the quantum system of interest is not coupled 
to any bath but coupled to an IHO which is coupled to a bath. The Hamiltonian of this 
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model reads 



H b = $r + U (q) + ^7 + "MO* (X + Ag) 2 



2/i ^' 2M 2 



In the model, the pure Hamiltonian of the quantum system of interest is p 2 /2fi + U (q). The 
IHO has mass M and frequency f2 , and is coupled by the quantum system of interest with 
strength A. The bath coupled to the IHO is comprised of a set of harmonic oscillators, and 
their mass, momentum, coordinate, and coupling coefficients are denoted by {mj,pj, x i: q}. 
The introduced Ki is a controlling parameter for discussions. It has been shown that the 



model B can be reduced to the model A through an ESDF of the bath 
cut-off frequency of the bath be infinite, we can obtain the ESDF as 



7]. Setting the 



j| M = , (17) 



which has been obtained in Ref. Q. Here, T = Kif]/M. If the cut-off frequency is finite, the 
ESDF becomes 

j|M= A ^ MM 2 . (18) 

Here, 

Z(cu) = -u 2 + tt 2 + KiuTRelW (to)]. (19) 

If the quantum system of interest is a qubit the model B can be reduced to a spin-boson 
model as Eq. (fl5i) with ESDF Jb {uj). By using the ITM algorithm we can also investigate 
the reduced dynamics of the qubit. The evolutions of the diagonal and off-diagonal elements 
of the reduced density matrix of the qubit governed by J B {uj) and (cj) are plotted in 
Fig.3 for the cases of u c = 3. OA, 5. OA, 10. OA, and 25. OA. Here, we set Q = 52. OA, and 
r/ = 0.0035, and other parameters are the same as Fig.2. It is contrary to the case of model 
A, here the qubit has longer decoherence and relaxation times governed by the J B {oo) than 
by Jf (u). 

As the cut-off frequency of the bath is not very large, the J\ (uj) results in shorter deco- 
herence and relaxation times for the qubit in the models A than the (uj) does. However, 
in model B, to compare to the J B (uj), the J B (u) causes longer decoherence and relaxation 
times of the qubit, as the cut-off frequency of the bath is not very large. When the cut-off 
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frequency of the bath increases to 5 times (for the model A) and 25 times (for the model 
B) of the characteristic frequency of the qubit, the J F /b (w) and the J^ B (ui) will lead to 
almost the same reduced dynamics of the qubit for the respective models. 



B. Two developed models (models C and D) 

In this section we investigate the other two models C and D. These two models are 
developed from above two basic models A and B. 
Model C: The Hamiltonian of the model C reads 

2 1 / v\ 2 

Pi , 1 2 / , KlCiX X 



2m,- 2 V rriiUJ 



+ ^ + M^ + 5f) 2 ' (20) 

Here, the quantum system of interest interacts with the IHO and both of them are coupled 
to their independent baths. The bath {B\) coupled to the IHO and the bath (B2) coupled to 
the quantum system of interest are considered to be constructed with two sets of harmonic 
oscillators and those mass, momentums, coordinates, and coupling coefficients are denoted 
by {m a , p a , x a , c a } (a = i for Bi, and a = j for B 2 ). As k±, the k 2 is also a controlling 
parameter for discussion. Similar to above subsection we can obtain an ESDF for the model 
as 

J 1 H = V UK 2 2 e-^ + 2 X2Q o^V (21) 

when the cut-off frequency of the bath is infinite. Setting the cut-off frequency of the bath 
be finite, we have the ESDF for the model as 

J F C (u;) = tjukIQ (u) + — -4— \ ; \ (22) 

If both A and K\ are equal to zero, namely the quantum system of interest have not any 
interaction to other systems except for being directly coupled to a bath (where the bath is 
B 2 ) the model C will reduce to the model A, and from Eqs.(21) and (22), we can obtain 

jy F (u) ^> JT H ■ (23) 
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When k 2 = 0, the model C reduces to the model B which describes a quantum system of 
interest interacting with its environment through an IHO, and from Eqs.(21) and (22), we 
have 

jf( W )^jf( W ). (24) 

Model D: Model C describes that the quantum system of interest interacts with an IHO and 
both of them are coupled to two independent baths. However, when the quantum system is 
very closed to the IHO, they may embed in one common bath, then the Hamiltonian of the 
total system becomes 



2 9 

,2 



2 

Pi 1 2 



Xi H (kiX + K 2 q) 



rriiUjf 



(25) 



l 2m; 

Similarly, we can obtain an ESDF from infinite cut-off frequency as 

J* D (w) = \MQ 2 Im $ + cor] [ Kl K2 Re $ + 4] e - w/Wc . (26) 

Setting the cut-off frequency of the bath be finite, we have 

j£ ( w ) = \Mtt 2 Im + urjK^ Re W (u) Im 

+ wr? [/cj/caRetf + k|] 8(u), (27) 

where 

(cu 2 - Q 2 ) + m^Ve-^ ' 1 ; 

If we set Ki — 0, A = 0, and k 2 — 1, the model D will reduce to the model A and from 
Eqs.(26) and (27), we have 

T I/F I \ A=K1=0,K2=1 T J/F / \ / Qn x 

M > J A ( w ) • ( 3 °) 

Setting k 2 = 0, the model D will reduce to the model B and from Eqs.(26) and (27), we can 
obtain 

J^H^J^H. (31) 
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By using the ESDFs f c ft (to) and J]** (u) we can compare the reduced dynamics of the 
qubit in models C and D, as we have done in last subsection. We can plot above different 
spectral density functions as Fig. 4. Here, we are more interested in the question: in which 
model (C or D) the qubit has longer decoherence and relaxation times when the values of all 
of the parameters in the models are the same? Using the spectral density functions J c (uj) 
of the model C, and J^ F (co) of the model D, we solve the reduced dynamics of the qubit 
in the two models. The evolutions of the elements of the reduced density matrix of the 
qubit are plotted in Fig. 5. Here, we set Ki — k,2 = A = 1, fi — 10A, and uo c = 7A. Other 
parameters are set with the same values as in Fig.2. It is shown that the qubit in the model 
D has longer decoherence and relaxation times than it has in the model C at the same finite 
cut-off frequency of the baths. 



III. DISCUSSIONS AND CONCLUSIONS 



The spectral density function is in fact a bridge to connect the microcosmic and macro- 
scopical physical quantities in dissipative quantum theory. At low temperature, the exci- 
tations of the bath modes with high frequencies are suppressed. So the spectral density 
function derived from infinite cut-off frequency of the bath will be out-of-true. In this paper 

TIT? 

we investigate the ESDFs J x (oj) for models (x =A, B, C, D) at infinite and finite cut-off 
frequencies of baths from the microcosmic dynamical models by using the method initiated 
by Leggett et. al. As the cut-off frequency of the bath is comparable to the characteristic 
frequency of the quantum system of interest, the non-Markovian effects may also be impor- 
tant. In this case, it is improper to use any Markovian approximations to investigate the 
the reduced dynamics of quantum system. So we used a numerically exact ITM algorithm 
based on the QUAPI scheme solving the reduced dynamics of qubits. 

In order to clearly compare the difference of the dynamics of quantum systems resulting 
from different ESDFs (u) and (u), we used the effective coupling constants calculating 
the dynamics of the qubit in each models and fixed other parameters such as the tempurature 
is set 300k. It is shown that a qubit (the quantum system of interest) in these models has 
different dynamics when they are governed by the spectral density functions (u) and 
Jjf (oo) as the cut-off frequencies are not too large. The difference increase with the decrease 
of the cut-off frequency. When the cut-off frequencies of the baths are larger about 5 times 
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(for the models based on model A) and 25 times (for the models based on model B) than 
the characteristic frequency of the qubit, the two kinds of spectral density functions result 
in almost the same reduced dynamics. In addition, we obtain that a qubit will have longer 
relaxation and decoherence times when it is coupled to an IHO and both of them are coupled 
to a common bath rather than to their respective baths. 

Acknowledgement 1 This project was sponsored by National Natural Science Founda- 
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No. 2008 A6 10098) and K C Wong Magna Foundation in Ningbo University. 



IV. APPENDIX 

In the Appendix we give a proof of Eq. ffTTj) . namely prove equation R (cu) = -n e -^/ M _ 
l~W (u) . Setting 

z = x + i%f(z) = ^M, (32) 

(z z — cu z ) 



and from Fig.5 we see that there is not any singular point in the loop of (L + M2O + OM1 +C). 
So by using the Cauchy's theorem we have 

f(z)dz+ f f{iy)d{iy)+ I* f{x)d{x)+ [ f(z)dz = 0. (33) 
Jr Jo y Jc 

Gi G11 Gin=R(u>) Giv 

For curve L, because 

A x = lim zf(z) = hm ^tlM = (34) 

we have 

d = J f(z)dz = i (fc - 0i ) Ai = 0. (35) 

For curve C, we replace pramaters z and u with r and cp and setting z — u = re lip . Because 
ip 2 = 7T, <pi = 0, and 

A 2 = lnn z/W = ^M, (36) 



we have 



Giv = I f{z)dz = % {(p 2 - fi) A 2 = — . (37) 
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For linear M2O, the integral becomes 



Setting - = x, — = m, we have 



Here, 



i f° cos rra , . 1 sin mx 
= / d(z) / d (z) . 



cos mx 

^ COS 7TIX 7T 

a (x) = — — exp(— 0J/uj c ), 



x 2 + l y ' 2 



sin mx , . . 
a ix) 



R 



with 



x 2 + 1 

— Shi (m) cosh (m) + Ci (— im) sinh (m) + — sinh (m) 
ttH/ (a;). 



. sinh(t) , 
Shi{m) = / —dt, 







Ci(m) = /°° e'H^dt + ln(m) + -d*. 



So we have 



G n = — e^ c + -W(u) 
2u ui 



From Eqs.f l33p . (I35I) . (137|) . and (JHJ) we have the integral in linear OMi as 



Gjjj= / /(x)d(x) 



— G/ — Grj/ — Gjy 

2uo 2tu 2u y ' 

_V wc --iy( w ) 

U CO 

R(u). 
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FIG. 1: (Color online) The values of Im[VF (w)] when choosing the cut-off frequency oj c different 
values. It is shown that the absolute values of Im[lF (w)] decreases with the increase of the cut-off 
frequency of the bath uj c 
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FIG. 2: (Color online) The reduced dynamics of qubit governed by j\ (w) vs by j\ (oj) in different 
values of ui c . Top four figures a, b, c, and d plot the evolutions of the diagonal element pn{t) 
beginning at the initial state /5i(0), where A = 10 12 Hz, e = 0.01 A. Bottom four figures e, f, g, 
and h plot the evolutions of off-diagonal element pi2{t) beginning at the initial state P2(0), where 
e = 10 12 Hz, A = O.Ole, T=300 k, = 0.004. The pi(0) and p 2 (0) are supposed in the text. Here, 
and in the following, we denote the results from the ESDFs J 1 (co), and J F \lo) with I and F in the 
Legend 
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FIG. 3: (Color online) The reduced dynamics of qubit in model B governed by J l B (oj) vs by Jg (oj) 
in different values of oj c . Here, A = 1, K\ = 1, rj = 0.0035, and T = 52A, and the values of other 
parameters are the same as in Fig. 2. 
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FIG. 4: Effective spectral density functions, J% (oj) in different models x = A (black asterisk), 
x = B (blue diamond), x = C (red circle), and x = D (green box), and J* (oj) in different models 
x = A (black cross), x = B (blue solid diamond), x = C (red solid circle), and x = D (green solid 
box). Here, we set r] = 0.02, A = 1, K\ = k 2 = 1, T = 52A,f2 = 10A,a; c = 11A,cj = 10A 
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FIG. 5: (Color online) The reduced dynamics of the qubit in model C versus in model D governed 
by the ESDFs J c (oj) and (w). Top two figures a and b plot the evolutions of the off- 
diagonal pi2(t)- Bottom two figures plot the evolutions of the diagonal elements pn(t). Here, 
K\ = Ki = X = 1, r/ = 0.0035, and u c = 7 A, S7o = 10A, and the values of the other parameters are 
the same as in Fig 2. 




FIG. 6: Integral scheme of Eql33l 
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